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Combining the self-consistent theory of localization and the dynamical 
mean-field theory, we present a theoretical approach capable of describing 
both self-trapping of charge carriers during the process of polaron forma- 
tion and disorder-induced Anderson localization. By constructing random 
samples for the local density of states (LDOS) we analyze the distribution 
function for this quantity and demonstrate that the typical rather than 
the mean LDOS is a natural measure to distinguish between itinerant and 
localized states. Significant polaron effects on the mobility edge are found. 
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The question of how the electron-phonon (EP) interaction influences 
the localization transition caused by disorder [1], i.e. by strong impurity- 
induced spatial fluctuations in the potential energy, has been addressed by 
Anderson about thirty years ago [2]. He called attention to the particular 
importance of EP coupling effects in the vicinity of the so-called "mobility 
edge", separating itinerant (extended) and localized states. Nevertheless, 
there is as yet not much theoretical work even for the simplest case of a 
single electron moving in a disordered, deformable medium. 

As a first step towards addressing this problem, in Ref. [3] the single- 
particle Holstein model with site-diagonal, binary-alloy-type disorder was 
studied within the dynamical mean field approximation (DMFA) [4]. The 
DMFA, however, cannot (fully) discriminate between itinerant and localized 
states, mainly because the randomness is treated at the level of the coherent 
potential approximation. In order to remedy this shortcoming, recently the 
authors [5] adopted the statistical DMFA (statDMFA) [6] to the Anderson- 
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Holstein Hamiltonian, 

n = ^ - J 5>k- + H - c -) - \/^X>* + b b n i + n E b ^ ' 

i (ij) i i 

where J denotes the electron transfer amplitude, f2 is the frequency of 
the optical phonon, E p is the polaron shift, and the on-site energies {e^} 
are assumed to be independent random variables with probability density 
p(ej) = (1/7)^(7/2 — |ej|). The statDMFA is essentially a probabilistic 
method (in the sense of the self-consistent theory of localization [7] ) , based 
on the construction of random samples for the physical quantities of interest. 

As a natural measure of the itinerancy of a polaron state, we consider 
the tunneling rate from a given site, defined - on a Bethe lattice with connec- 
tivity K ( J = J\[K~) - as the imaginary part of the hybridization function 

K 

T 1 {uj) = (ttJ 2 /K)Y,Ni{lo), where N^u) = -(l/^ImftH (2) 
1=1 

is the local density of states (LDOS). The LDOS, directly connected to 
the local amplitude of the electron wave function, undergoes a qualitative 
change upon localization implying a vanishing tunneling rate Ti(uj) for a 
localized state at energy to. The local single-particle Green function and 
the related hybridization function are given by (z = 00 + irf) 

1 J 2 K+1 1 

Gi{z) = z-ei-HM-Xiiz) ^ Hi= K^ [ z-e l -Hi-^ 

(3) 

respectively. We now ignore that the functions on the rhs of Hi should be 
calculated for the Bethe lattice with the site i removed, i.e. we make the 
replacement {Q\, H\, T,]} ~> {Qi, Hi, £/}, and furthermore take K as 
the typical number of terms even for the central site. Finally, the EP self- 
energy contribution is determined in the limit K — > 00. The self-energy is 
then local and, in terms of a continuous fraction expansion, takes the form 

„ / n E n 1 , , . 

{Z)l [if>(*)]-i-... 

with [fI p) (z)]- 1 = z- P n-ei-Hl P \z) and H { f> \z) = H^z-pO). Here the 
energy shift keeps track of the number of virtual phonons (0<p<M). Re- 
gardless of the local EP self-energy, the statDMFA takes spatial fluctuations 
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of, e.g., the LDOS into account and provides an adequate description of dis- 
order effects. Due to the randomness in the on-site energies, the tunneling 
rate and consequently the LDOS is a random variable, and the question 
of whether it vanishes or not depends on the probability density exhibiting 
different features for itinerant and localized states [1, 7]. In particular, the 
difference between the mean and typical LDOS, 



obtained by the arithmetic and geometric mean of the LDOS, respectively, 
is a useful measure to discriminate between extended and localized states. 
iV mean (w) > but N typ (uj) = indicates a localized state at energy oj. 

In the numerical work, we calculated the LDOS by solving a recursion 
scheme for which depends on Ke's, KH^ 's, . . . , and KH^ Pmax) 's. 

Starting from an initial random configuration for the independent vari- 

(p) 

ables , which is successively updated with a sampling technique similar 
to the one described in Ref. [7] , we constructed self-consistent random sam- 
ples for h[ p \ using K = 2, N = 100 000, M = 35, and rj = 10~ 8 . 

Without disorder, the physical properties of the Holstein model are de- 
termined by two interaction parameters, A = E P /2J and g 2 = E p /Q, and 
the adiabaticity ratio a = VI/ J. Polaron formation sets in provided that 
A > \/\[~K and g 2 > 1. Of course, the internal structure of the polaron 
depends on a. Disorder affects polaron states quite differently in the adi- 
abatic (a <C 1), non-adiabatic (a ~ 1), and antiadiabatic (a S> 1) cases. 
Without EP coupling, i.e. in the pure Anderson model, the critical disor- 
der strength needed to localize all states is (7 c /W / o) complete ~ 2.25, where 
Wq = 4J. In the weak EP coupling regime, it has been shown that the 
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Fig. 1. Mean and typical LDOS for 
the lowest polaron subband in the 
non-adiabatic strong EP coupling 
region (A = 9.0, a = 2.25, and 
J = 0.5). The pronounced disorder- 
induced broadening of the LDOS 
occurs because the variation of the 
on-site energies 7 = 0.001 is on the 
order of the strongly renormalizcd 
band width W. 
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Fig. 2. Left panel: Mean and typical LDOS in the adiabatic intermediate-to-strong 
EP coupling region (A = 1.0, a = 0.25, J = 0.5). Note that N mcan (u) is almost 
perfectly approximated by the DMFA. At about u) = —1.13 the second polaronic 
subband starts. Right panel: Probability density of the LDOS for four representa- 
tive energies uj. The inset shows the probability distribution, i.e., the cumulant of 
the probability density. 



quantum interference needed for localization is significantly suppressed by 
inelastic polaron-phonon scattering processes [5]: States above the optical 
phonon emission threshold are more difficult to localize than the correspond- 
ing bare electron states. In the very strong EP coupling regime, extremely 
weak disorder turns itinerant into localized polaron states. Surprisingly, the 
ratio (7c/W)compiete, where W is the band width of the lowest polaron sub- 
band, is almost the same as for a bare electron. In fact, in the non-adiabatic 
strong EP coupling regime, where the band collapse changes only the overall 
energy scale, disorder affects a polaron in a similar way as a bare electron. 
For example, the LDOS and mobility edges are symmetric (cf. Fig. 1). 

In the adiabatic intermediate-to-strong EP coupling regime the physics 
is much more involved. Here the band dispersion of the lowest subband 
significantly deviates from a rescaled bare band [8], leading to a strong 
asymmetric LDOS. Specifically, the states at the bottom of the subband 
are mostly electronic and rather mobile due to long-range tunneling induced 
by EP coupling, whereas the states at the top of the subband are rather 
phononic and immobile [8]. As a direct consequence, the states at the 
zone boundary are very susceptible to disorder, i.e. the critical disorder 
strength needed to localize these states is much smaller than for states at 
the bottom, and, from the results for the typical LDOS, we find asymmetric 
mobility edges (see Fig. 2). Moreover, (7 c /W) CO mpiete ~ 2.8, which is larger 
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than the corresponding ratio for a bare electron. Thus, contrary to naive 
expectations, at intermediate EP couplings, an adiabatic polaron is even 
more difficult to localize than a bare electron. 

It is very instructive to discuss the behaviour of the probability density 
of the LDOS and the corresponding probability distribution. Note that 
both quantities have to be calculated self-consistently within our sampling 
procedure. The right panel of Fig. 2 shows the dramatic change of the 
probability density of N[(uj) when the system undergoes the localization 
transition by crossing the mobility edge. In the region of localized states, 
the probability density for the LDOS is broad and very asymmetric and, as 
a consequence, the mean LDOS is not representative. 

In conclusion, in terms of the Anderson Holstein model, we have demon- 
strated that the statDMFA, which according to the spirit of Anderson's early 
work [1] focuses on distribution functions and associates typical rather than 
mean values to physical quantities, yields a proper description of disordered 
electron-phonon systems. 
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